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Abstract 

After a short exposition of field correlators in the QCD vacuum 
and the recently discovered Casimir scaling phenomenon, the origin 
of confinement in QCD is discussed and two possible mechanisms are 
suggested, which can be checked by new lattice measurements. Screen¬ 
ing of confinement due to sea quarks is discussed and quantitatively 
explained. Deconfinement is introduced via the colorelectric field evap¬ 
oration and the transition temperature Tc is found numerically in good 
agreement with lattice measurements. The dependence on and 
n/ is also predicted and agrees with the recent lattice data. 


1 Introduction 

The confinement is known to be the most important QCD dynamics at large 
distances preserving stability of matter and existence of our world (for a 
review see [Q). The force of approximately 15 tons between quark and an¬ 
tiquark in mesons, or between quark and the string junction in proton is 
known from the lattice with excellent accuracy to be constant up to very 
small distances P|. However theoretical understanding of this phenomenon 
still far from complete despite many efforts during the last decade. It is a 
purpose of this talk to describe in some detail what is understood now in 
the picture of conhnement and deconhnement and where is the front line of 
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our present knowledge. The main instrument in what follows is the gauge- 
invariant Formalism of Field Correlators (FFC) P, P], however references to 
particular applications are done at appropriate places. 

As it is understood now the nonperturbative structure of the QCD vac¬ 
uum responsible for mass generation and conhnement is adequately described 
by gauge-invariant held correlators 




where <F(a;, y) = Pexpig A^dz^ is the parallel transporter, and (...) means 
the vacuum average with the standard QCD action. 

The ensemble of D^'^\ n = 2, ...oo, = 0) or equivalently, the ensemble 

of connected correlators (called cumulants) contains a complete infor¬ 
mation for the quark-antiquark dynamics in the limit of large Nc (at hnite 
Nc also another set of correlators containing more color traces is necessary). 

A recent lattice data p for Wilson loops in different representations of 
SU(3) color group have conhrmed the notion of Casimir scaling with very 
high accuracy (around 1%), i.e. the fact that static QQ potentials are pro¬ 
portional to the quadratic Casimir operator. As it was argued in p it is 
the lowest held correlator which has the property of Casimir scaling, 
while all higher correlators violate this property and hence are strongly sup¬ 
pressed in the vacuum. This remarkable observation makes the picture of 
the QCD vacuum rather simple, it can be called the Gaussian Stochastic 
Vacuum (CSV) and it also enables one to use the lowest correlator known 
from lattice calculations 0] for all nonperturbative dynamical calculations, 
including the study of QCD strings in mesons P and baryons 0, study of 
spectra of mesons, hybrids, glueballs and baryons (for a review and references 
see 0, 0). 


At the same time one should stress that to make the theoretical approach 
of FFC complete, one should calculate FC in the framework of the same 
method using the only QCD parameter, (string tension a or Aqcd), as an 
input, and in addition to understand the dynamics of Casimir scaling, i.e. 
suppression of higher cumulants. 

In this direction only hrst steps are done 


with encouraging results. 


in particular the smallness of the gluonic correlation length Tg, entering FC, 
is understood from the gluelump spectrum [|T^]. 

However, the understanding of the Casimir scaling is not yet complete. 
On one hand, the contribution of higher correlators to the string tension 
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in general can be estimated as {FTg)^~^F, where F is average field strength, 
estimated from the condensate to be F ~ 0.5 GeV^, and here we have not 
taken into account that cumulants are connected. If this is included, then 
the Casimir scaling violation can be shown to arise due to white exchanges 
between ’’dipoles” {F^F), while higher cumulants are suppressed as 1/N^ 
(see last ref. in P). Altogether one gets a suppression factor ~ 

3% which gives an order-of-magnitude agreement with the 
violation of Casimir scaling on the lattice 0. 

Another set of relevant questions concerns the nature of confinement, i.e. 
which field configurations are responsible for confinement. In the GSV it is 
the lowest correlator which confines, therefore the question is about the 
nature of field configurations which saturate D^‘^\ It is the purpose of this 
talk to discuss these points one by one, starting with confinement mechanism 
in the next section, and temperature deconfinement in section 3. 


2 Mechanism of confinement 


Since according to lattice data ensures some 99% of confinement, one 
should look more carefully at its structure, @], namely 






+ {fii ^ Ui)]Di{z), z = xi-X2, df,= 


Zu2^Ui^2)F 

d 


dz„ 


( 2 ) 


Using the nonabelian Stokes theorem and the cluster expansion (see 
for details and references) one has for a large Wilson loop 


{W{C)) = ^(trPexpi^i f A^dx^) = eyzp{-aSmin) (3) 

aVc Jc 

where 

~ 2 J F{u)d^u (4) 

and Smin is the minimal area inside the contour C. From (^ it is clear 
that D{z) plays the role of the order parameter for confinement (at least for 
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Nc —>■ oo) and this is confirmed by lattice calcnlations of D{z), where D{z) 
vanishes abrnptly above the critical temperatnre T,, 0]. A farther analysis 


can be done applying to both sides of 
yields [|T4 


the operator which 


dD{z 


dz‘^ 


1 = f(b 4- f(2) f(b = 
J a 'Jo.') J a — 


24:N, 


[{trDyF^p{z)^{z, 0)Fap{0)^{0, z))]] 


2 pz 

J ‘iysO‘(y){trFxp{z)ii{z,y)Fxs{y)t(y,0)Fcfj(0)ti(0,z))-h.c.] 

( 5 ) 

We start with the Abelian case where one shonld replace d^, 1 

in d^) and the last two terms inside sqnare brackets, coming from the contonr 
differentiation, are absent: 


Qj^ahelian 

dz"^ 


= /t’ = ainst{i']l,'“"‘\z)F„i,((l)) 


( 6 ) 


and = d.yFzfj 3 [z) is the monopole cnrrent. Hence in the Abelian U(l) 

case conhnement might be dne to Abelian monopoles and this fact agrees 
with what is known abont the Abelian mechanism of conhnement on the 
lattice [|I|. 

In the nonabelian case, however, the corresponding Bianchi Identities 
(BI) 

= Ju{x), (7) 


where J,y{x) = 0, are usnally assnmed to hold. In the Abelian projection 
method (see [|^ for a discnssion ) one separates from F^^ in some special 
gange a singnlar term, the color-diagonal monopole-type conhguration which 
violates both Abelian and nonabelian BI. As it is clear from d), (0 this 
procednre clearly snpplies a sonrce fa^ which makes D{z) nonzero. One 
shonld, however, make two remarks in this connection. Firstly, when one 
tries to associate the conhning conhgnration with some classical monopole, 
there appear no magnetic monopoles with hnite selfenergy and stable |T^ . 
Therefore in it was conclnded that only qnantnm conhning conhgnration 


can snrvive, which however not described analytically. 

Secondly, and this is a more fnndamental difficulty, associating conhning 
conhguration with necessary violation of BI, one usually abandons the con¬ 
nection of Ffj^tz with Afj_, F^y{A) = d^Ay — dyA^ — ig[A^,Ay] adding to this 
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expression extra terms (see e.g. in |^]). In this way one finds confinement 
not for the original QCD Lagrangian, bnt for a modified form. 

In what follows we shall insist on keeping the original expression F^y{A) 
intact. In this case the nonabelian BI reduces to the equation 


const 0 . 


( 8 ) 


One can see that violation of BI, Eq. (^, requires a very special form of 
Ap{x). One possible form of Ap[x) is given by the path ordered exponents, 

e.g. 


^ 0 {x) = Pexp[ 


<c 


dz^X^{z)]h 0 {x) + h.c. 


(9) 


with noncommuting fi v and the contour O, which belongs to a wide 

class of nondifferentiable contours; x is the end point of the contour, and P 
orders Xfj,{z) from say, right to left with 2 ; approaching x. Such configurations, 
if possible, would make nonzero Jy{x) and contribute to the first term on the 
r.h.s. of d^), thus yielding nonzero D{z) and string tension. In this way one 
would find the first possible source of confinement ~ due to the term fa^ in 
(|^). If this mechanism is proven, it would physically imply a monopole-like 
mechanism of confinement. One should stress that this does not mean real 
magnetic monopoles or magnetic fluxes present in the vacuum, since both 
violate Casimir scaling p, while D{z) supports it. The second source is due 
to the term fa^ in (^, which is always present irrespectively of the violation 
of the BI. To make explicit the meaning of this term one can take the limit 
z —>■ 0 and obtains Hi 


dD{x'^ 

dx"^ 
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r“{F“^(0)F»,(0)F'„(0)). 


( 10 ) 


Since the triple correlator can be written as (Ef(0)Ej'(0)i?^(0)) 

one can view this second mechanism as the creation of magnetic fluxes from 
electric fluxes, i.e. the electric fluxes contained in the parallel transporter 
when shifted in the process of differentiation, serve as a source of magnetic 
flux, thus replacing magnetic monopoles. 

To distinguish between these two possibilities one could measure on the 
lattice the correlator 0 


= j^tr{Df,Fi,^{x)^{x,y)DaF^p{y)^{y,x)). 


( 11 ) 
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The nonzero answer for would mean that the violation of BI is indeed 
the source of confining configurations and writing 

= {dadfj - d^6o,p)D{x) + aJJ], (12) 

where is the contribution of fa \ one can estimate the role of the first 
mechanism as compared to the second. Here Aap is measured via (pT]) and 
D{x) is measured directly from (^). 

One should note that the given above formulation of confinement mech¬ 
anisms is fully gauge invariant and does not need gauge fixing. 

It is worth mentioning that the FFC and the applicability of the method 
are not directly related to the question of the mechanism of confinement, 
since FFC exploits field correlators as input, but in the analytic calculation 
of D{x) and Di{x) the problem of confining configurations becomes essential. 

3 Deconfinement 

Discussion of this phenomena can be done in 3 different directions: 1) when 
Nc = 3 and not infinite, the string connecting static quark antiquark breaks 
up at some distance Rb (typically Rb ~ 1.4 fm) forming two heavy-light 
mesons. The same happens for a light qq pair; 2) when the temperature T 
exceeds Tc the color electric string between Q and Q disappears; 3) when 
baryon density pB exceeds critical value pc, confinement is believed to disap¬ 
pear. We shall discuss below only points 1) and 2). 

1) For finite W, e.g. W = 3, the effect of sea quarks is given by the 
quark determinant det(m -|- D) present in the integral measure of averaging 
the Wilson loop in (j^). This determinant is producing additional light-quark 
loops due to the heat-kernel representation 

1 /•OO 7 

(det(m-F£))fF(C')) = (exp[--fr / —{Dz)^^e~^w„{C^^)]W{C)). (13) 

z Jo s 

Expanding the first exponential in (|^ one obtains corrections to the static 
QQ potential due to the sea-quark loops WaiCxx) where Cxx is the closed 
contour passing through the point x, which is integrated over in (|T^) . 

The leading correction is proportional to the loop-loop correlator 

x{Ci,C2) = {W{Ci)W{C2)) - (W(C'i))((W(C'2)), (14) 
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which was calculated in [IS , and the energetically lowest conhguration cor¬ 
responds to the creation of a hole in the largest quark loop. These holes 
produce a partial screening of the QQ static potential and consequently a 
signihcant decrease of highly excited meson masses. A detailed analysis of 
this situation was done in where the form of the screened potential was 
found yield a perfect agreement of predicted light meson masses with exper¬ 
iment (for L = 0 ,15 2, 3 and = 0,1, 2, 3,4). 

Note that this mechanism of the screened conhnement is different from 
the two-channel model, since in the former sea-quark loops are virtual and 
do not cause actual decay. Physically this difference results in the relative 
insensitivity of the meson mass depletion on the channel quantum numbers 
and on the proximity to the decay threshold in the screened conhnement 
picture in contrast to the two-channel approach. 

2) The temperature deconhnement 

We start with gluodynamics and use background perturbation theory for 
T > 0 1201 which allows to separate in = B 


^ background gluon 
held from valence gluons a^, both with periodic boundary conditions at 
X 4 = n(3 = n/T. 

The main ehect of nonzero T is that the correlators of electric helds, 
g‘^{Ei{z)Ek{0)) = 5ikD^{z) + 0{Df), and magnetic helds g'^{Hi{z)Hk{0)) = 
6 ikD^{z) + 0{D^) are diherent, D^{z) 7 ^ D^{z), Df{z) 7 ^ D^{z). The same 
is true for the condensates, (E^) 7 ^ (H^), which dehne the vacuum energy 
density 

£0 =-llw—(E' + H^). (15) 

“ 3 167r' ^ ^ ^ 

It is clear that the conhnement phase of gluodynamics with e = ^(l) 
consists of glueballs moving in the conhning background (E^ + H^), while 
the deconhnement phase with e = e{ 2 ) represents valence gluons moving in 
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the deconhned vacuum background (H"^). It was understood already in 
that it is advantageous (for the minimum of the Free Energy of the Vacuum 
(FEV)) to keep magnetic condensate intact at T > Tc, while at T < Tc both 
electric and magnetic can be nonzero. This implies that and the spacial 
string tension agp stay constant across phase transition which was supported 
by lattice data [0]. One can dehne FEV in two phases including quarks as 

0 

TT 

(higher mesons) 


TT 


F(1)/U = £(1)-^T* 
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F(2)IV, = e(2) - (Nl - 1 ) 




TT 


45 




Ttt^ 

im 


NcT^UfVtq + 0{N[ 


(16) 


where Hg, VLg are perimeter contributions to the quark and gluon loops respec¬ 


tively. Assuming that magnetic part of £(1) does not change for 0 < T < 
and equal to the electric one, one has £(1) = Sq. e{2) = ^Eq and from 
F(l) = F{2) one obtains for the transition temperature [p^l : 


T. = 


_Ifol_ 

;=g-(NI - m, + > 2 ), 


1/4 


(17) 


Neglecting the influence of magnetic helds on the quark and gluon motion 
yields fig = fig = 1, and one can compute for different Uf taking Eq 


([T5| ) from standard gluon condensate estimate 
Uf = 2, 3,4 and 3 times more for ny = 0 
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= 0.012 GeV^ for 
One obtains in this way from 
([TTp, Tc = 240,150,141,134 MeV for rif = 0, 2, 3,4 which can be compared 
to the lattice data = 270,172,154,131 MeV respectively. One can see 
a systematic 10-12% disagreement which can be removed taking standard 
gluonic condensate 1.5 times larger. 

Note that ([T7|) gives for Tc the constant values for large Nc, and Tc is only 
weakly dependent on Nc- 

Moreover, the phase transition is the hrst order for large N^. Both facts 
agree with recent lattice data 
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Indeed the analysis of lattice data at 
^0, and for a = 0.18 GeV^ it gives = 


n/ = 0 in yields: 0 = 0.582 = 

(0.246 -F 0.02^) GeV, which agrees well with (|3) for fig = fig = 1, where 
Tc = 0.243 GeV and forn/ = 0 it does not depend on Nc, so that (signihcant) 
Nc dependence appears only for Uf > 0. The simple picture described above 
is sometimes called ’’the Vacuum Evaporation Model” (VEM) and is actually 
the only known model successful in all these features. However, the latent 
heat for Nc = 3, nj = 2 comes out in VEM too high, and to improve the 
situation one should take into account that fig, fig ^ 1 above Tc and the 


contribution of other than pion mesons in (|T^) below Tc, which was done 
Both factors make the phase transition more smooth and strongly 
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decrease the specihc heat. 

In conclusion, it is clear that the present approach yields a very consistent 
picture of both conhnement and deconhnement, which can be derived from 
the standard QGD Lagrangian using background perturbation theory and 
FFG. In this approach inputs are D{x), Di{x), which are known from the lat¬ 
tice data 0 to be D{x) = Z1(0) exp(—|a;|/Tg), Di{x) = Zli(O) exp(—|a;|/Tg) 

































reducing input to 3 numbers, Zi)(0), Zi)i(O) and Tg. Since -Di(O) -D(O), one 
actually has -D(O) and Tg, or equivalently a and Tg, and most hadron spec¬ 
trum data are computed through only a with 10% accuracy. Nevertheless 
the method can be considered as logically consistent, with the statement that 
everything including all held correlators is computed through the only QCD 
parameter: Aqqd or o'- The work in this direction was partly done in |^|, 
There the correlator length Tg was expressed through the gluelump 


masses and Tg was found in the range 0.13-0.17 fm. 

Moreover the full correlator D{x), Di{x) can be computed through the 
gluelump Green’s function [^, if one assumes that the nonabelian BI are 
not violated. In the opposite case there appears the problem of identihcation 
and study of the Bl-violating conhgurations, similar to the study of mag¬ 
netic monopole ensemble, if they really exist in the QCD vacuum. Indeed 
the Casimir scaling phenomenon strongly limits the admixture of any 
coherent conhgurations (like classical solutions for monopoles, dyons and in- 
stantons) with size larger than 0.2 fm Q, which makes the popular picture of 
the QCD vacuum with magnetic monopoles or central huxes rather unrealis¬ 
tic. Therefore for conhnement one has either the picture of the BI violating 
conhgurations with very small radius or an alternative picture with vanishing 
BI and the gluelumps saturating the held correlators. 

Present lattice data cannot unfortunately distinguish between two alter¬ 
natives since, as argued above, the popular abelian projection method or 
central vortex projection both measure the quantities proportional to D{x) 
or its derivatives, while only direct measurement of the correlator ([IID 
can separate two possibilities, provided a clean lattice analog of is for¬ 
mulated. 
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